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Abstract 

The reaction center of Chloroflexus aurantiacus and Rhodobacter sphaeroides mutation of RC's 
was investigated. To describe the kinetic of the Chloroflexus aurantiacus RC's we use incoherent 
model of electron transfer. It was shown that the asymmetry in electronic coupling must be 
included to explain the experiments. For the description of Rhodobacter sphaeroides H(M182)L 
mutation we used partially coherent as incoherent models of electron transfer. These two models 
are discussed with regard to the observed electron transfer kinetics. It can be concluded that 
partially coherent model is more adequate. We predict some new electron pathways for describing 
the kinetic of RC's and some mutation. 
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I. INTRODUCTION 



Photosynthesis is a reaction in which light energy is converted into chemical energy. The 
primary process of photosynthesis is carried out by a pigment-protein complex embedded in 
the membrane, that is, RC. In photosynthetic purple bacteria, the cyclic electron transfer 
reaction is performed by RC and two other components: the cytochrome (Cyt) hc\ complex, 
and the soluble electron carrier protein. The photosynthetic reaction centers (RC) is a special 
pigment-protein complex, that functions as a photochemical trap. The precise details of the 
charge separations reactions and subsequent dark electron transport (ET) form the central 
question of the conversion of solar energy into the usable chemical energy of photosynthetic 
organism. The function of the reaction center is to convert solar energy into biochemical 
amenable energy. Therefore, we wish to understand which features of the reaction center 
are responsible for the rate constants of these reactions. 

Insight into the molecular organization of the RC has been derived, initially, from spec- 
troscopic studies and, subsequently, from the development and analysis of high-resolution 
crystal structures of several photosynthetic organisms. The first RC structurally resolved 
(3 A) was of the purple bacterial RC from Rhodopseudomonas viridis This was soon 
followed by the elucidation of several other purple bacterial structures. Good progress is 
also being made toward achieving two- and three-dimensional structures of photosystem II 
(PSII) crystals. It is surprising that the structures of all of the different RC's show a dimeric 
core with a pseudo-C*2 axis of symmetry. 

A remarkable aspect of the RC structures is the occurrence of two almost identical elec- 
tron acceptor pathways arranged along the Ci axis relative to the primary charge-separating 
dimer (bacterio) chlorophyll (Fig. 1). This finding posed a key question: Does electron trans- 
fer involve both branches? In the purple bacterial RC, only one branch is active although 
the inactive branch can be forced into operation with modification of amino acid side chains 
on the active branch Moreover the charge-separating electron transfer reactions occur 
with a remarkably high quantum yield of 96%, where from two possible symmetric branches 
only the branch L is active in the electron transfer. This efficiency relies on the rates of 
the charge-separating reactions being 2-3 orders of magnitude faster than the rates of the 
competing reactions. 

The strong asymmetry imposed on primary charge separation photo-chemistry in the pur- 
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pie bacterial RC results from two homologous polypeptides that function as a heterodimer. 
A heterodimer is also involved in the core of the RC's of PSI and PSII. However, some RC's, 
such a S h e„ ob ac ter ia fl and gre en sultur bacteria fl, contain two identical —eric 
polypeptides, and electron transfer is potentially bifurcated. 

Genetic sequence information has greatly improved the understanding of the origin of 
the RC proteins. From the sequence analysis, it became clear that the purple bacteria RC 
is remarkably similar to that of PSII, and PSI was also discovered to have similarity with 
that of the green sulfur bacteria [5(. Recent structural comparisons between PSI and PSII, 
for example, show a distinct structural homology, which suggests that even these two RC's 
likely share a common ancestor [6j]. 

In purple bacteria, the electron rate is sensitive to the free energy difference between 
the excited state and the charge-separated state but not to the relative distribution of 
electrons over the two macrocycles of the donor. After extensive studies, the rate is now 
established to be critically coupled to the properties of the bacteriochlorophyll monomer 
that lies between the donor and bacteriopheophytin acceptor (Fig.l). The involvement of 
the bacteriochlorophyll monomer may give rise to multiple pathways for electron transfer [?| 
and can partially determine the asymmetry of the electron transfer along one branch [8]. We 
believe that the reason for asymmetric ET between prosthetic groups located on different 
polypeptides is a different molecular dynamics. Dynamics of atoms causes the change of the 
electrical potential fields and the conformational variations influence the mutual orientations 
between cofactors. Then the energy gap and overlap of electronic wave functions fluctuates 
as a result in the system. The net result is a different fluctuation of electronic energy 
levels on prosthetic groups and also a different fluctuations of the overlaps of the electronic 
wave functions on L and M branches. On the other hand the chain located on subunit M 
is inactive in ET and the highly asymmetric functionality, however, can be decreased by 
amino acid mutations or cofactor modification. We used this approach to explain the effect 
of individual amino acid mutation or cofactor modifications on the observed balance between 
the forward ET reaction on the L-side of the RC, the charge recombination processes, and 



ET to the M-side chromophores 
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13|. 



The theoretical models that describe the charge transfer in reaction centers using pa- 
rameters with clear physical interpretations. Some of these input parameters can not be 
deduced from independent experimental work. The information regarding the energetic pa- 
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rameters, the medium reorganization energies, the high frequency modes, and electronic 
coupling terms can be achieve with quantum mechanical computations. But until now these 
parameters which characterize the reaction centers are not available. And so we use the set 
of parameters which fit the experiments. Several sets of parameters were used to describe a 
charge transfer in the RC. A set of parameters based on molecular dynamics simulations [ijj] 
corresponds to a dominance of superexchange mechanisms for the primary ET reaction in 



RC's. Another set of parameters 
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16| was used to fit experimental data. This second 



set of parameters derives a dominant contribution from the sequential mechanism. The first 
set of parameters has the larger reorganization energies and the greater coupling factors. 
This set of parameters makes the ET rate much larger than it is found in the wild-type 
proteins. The possibility to find out the input parameters from theory is comparison of 
observed kinetics for different mutated reaction centers. The problem is that not always the 
impact of mutation on the input parameters is clear. In this paper we focus on the electron 
transfer in two RCs. First is the RCs of the green bacterium Chloroflexus aurantiacus. The 
second is the RCs of Rhodobacter sphaeroides H(M182)L mutation. It is believe that both 
purple bacterial RCs and RCs from C. aurantiacus have a similar structure 
in this work the set of parameters that characterized the observed L-side experimental ki- 
netics of wild-type (WT) RCs of Rb. sphaeroides very well. The Chloroflexus aurantiacus 
RC's and the RCs of Rhodobacter sphaeroides H(M182)L mutation have structural similar- 
ity but charge separation kinetics are different. Both these RCs contain BPheo pigment in 
M-branch in the position where BChl monomer is placed in the WT reaction center. In 
contrast with this structural similarity, the H(M182)L mutant reveal the electron transfer 
through the M branch, in the Chloroflexus aurantiacus RC's the M-branch is inactive. 



13] . We adapt 



II. THEORY 



We attempt to analyze the possibility that ET asymmetry can be described by model 
which assumes that there exists the vibrational modes of the medium which has a sufficient 
time for relax to the thermal equilibrium after each ET step. We start by considering an 
electron transfer system in which the electron has N accessible sites, embedded in a medium. 
We denote by \j) the state with electron localized at the jth site and j = 1, 2, N. The j 
and k sites are coupled by Vjk- The interaction of the solvent with the system depends on 
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the electronic states \j) by Hj. 



The total model Hamiltonian for the system and medium is 
H = H + V, (1) 



where 



N 



H = J2\3)[^- iT i+ H MjI 



(2) 



3=1 



N 



(3) 



j,k=i 



where Sj is the site energy The parameter h/2Tj has a meaning of the lifetime of the electron 
at site j in the limit of the zero coupling parameter. It can characterize the possibility of 
the electron escape from the system by another channel, for instance a nonradiative internal 
conversion or recombination process. 

The Hamiltonian describing the reservoir consisting of harmonic oscillators is 



Here, m a and oo a are frequency and the mass of the ath oscillator, and dj a is the equilibrium 
configuration of the ath oscillator when the system is in the electronic state The total 
density matrix p(t) of the ET system and the medium satisfies the Liouville equation, 





d tP (t) = - l -[Hp{t) - p(t)H^\ = -iLp(t). 



(5) 



In the interacting picture, 





The Liouville equation in the interacting picture has the following form: 



d tP (t) = -i[VKf)p/(f) -pi{t)Vi{t)\ = -iL(t) Pl (t), 



(7) 



where 




(8) 
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Here we denote the total trace, and the partial traces over the ET system and over the 
medium by Tr, Tr e , Tr^, respectively. By definition Tr = Tr®Tr e . The population on 
state \j) at time t is given by 



P j {t) = Tr{\j)<j\pQt)). 



(9) 



We assume that the vibrational relaxation is sufficiently rapid so that the system can relax 
to thermal equilibrium after each ET step. This assumption determines a choice of projector 
operator. The projector operator D acting on an arbitrary operator B in the Hilbert space 
of the total ET system and medium is defined by [18( 



(10) 



TV 

DB = J2Tr(\j)(j\B)p 3 \j)(j\, 

3=1 

where pj is the equilibrium medium density matrix in the state \j), i.e., 

exp(-Hj/k B T) 
Pi ~ TrQ expi-Hj/knT)' 

Using the standard projection operator techniques 19|,[20(] we can derive a generalized master 



(11) 



equation for the populations, 

2T 



n k=l J o 



h 

N rt 



+ J2 f W kj (t-T)P k (r)d 



r, j = l,...,N, j^k, 



(12) 



where 



W jk (t) 



\v Jk \ 2 



h 2 



Re< exp 



r,- + r A 



h 



exp 



i\ej - e k ) 
h 



e~^ at + n a e Wat - (2n a + 1)] 



(13) 



(14) 



Here, n a = [exp(huj a /kBT) — 1] 1 is a thermal population of the ath mode and 

Ejk = 2 maUJ a(dja — d ka ) 2 

is the reorganization energy of the ath mode when system transfer from state \j) to state 
\k). 



III. MODEL OF REACTION CENTER 



To describe the first step of electron transfer processes in the reaction centers we have 
used the 5-sites kinetic model of RC. 




FIG. 1: The RC of purple bacteria are composed of three protein subunits called L, M and H. 
Dimer P is describing by molecule 1. Cofactors in the subunits L are: 3 represent (BClil^) molecule 
5 (BPhj,) and 6 is (Ql) and identically in the subunits M (BCuIm) is describing by molecule 2 
and molecule 4 represent (BPh^). Cytochrom C serve as a source of electrons for reaction center. 

We designate the special pair P as site 1, the sites 2 and 3 represent the molecules BClil^ 
and BChl L , and the sites 4 and 5 then represent the molecules BPh M and BPh L (Fig. 1). 
We assume that we can neglect the backward electron transfer from quinone molecules and 
so we use the complex energies of 4,5 molecules of RC. Based on experimental observations 
of ET in RC, it is expected that bacteriochlorophyll play a crucial role in ET. In this 5-sites 
model we have assumed that ET in RC is sequential where P + BChl~ is a real chemical 
intermediate. The imaginary part of energy level 1 describes the probability of electron 
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deactivation to the ground state. We describe the ET in the Chloroflexus aurantiacus RC's 
by the following kinetic model 



or 

d t Pi(t) = ~( 1 ± + k 12 + k n )P 1 (t) 

+k 2l P 2 (t) + k 31 P 3 (t), 

d t P2(t) = -(A&i + k 2A )P 2 {t) + k 12 Pi(t) + k A2 P 4 {t), 

d t P 3 (t) = -{kaa + k^PzW + kuP^ + ksaPnit), 
2F 

dtP*(t) = -(-^+h 2 )P 4 (t)+k 24 P 2 (t), 
n 

or 

d t P 5 (t) = -(-^ + h 3 )P 5 (t) + k 35 P 3 (t). 



(15a) 
(15b) 
(15c) 
(15d) 

(15e) 



Here we denote fcy(s — > + ) = kij and kij(s — > + ) = J °° Wij(t)dt. 

We assume that the rate constant which characterizes ET can be described by both a low 
frequency medium vibrational mode and a high frequency intramolecular vibrational mode. 
At a high temperature regime the constant kij(s — > + ) is in the form 211 ] : 



h ■ - 



X 



/*°° 9tt 1 

( ^W* = T^ ( 4^-^f» I/2exp( - S ' 

ES ci j {Gji + A m ij + nhuj c ij) 2 , 



(16) 



Here, G^- = e.; — and S^y = ^m c ijUJ c ij(dci — d C j) 2 is the scaled reorganization constant 
for the high frequency ij-th mode, which is nonzero when electron is transferring from the 
state \i) to the state \j), and \ mi j = \m mi juj 2 n ^{d rn i — d m j) 2 is the reorganization energy of 
the low-frequency mode when the electron is transferring from the state \i) to the state 
The back electron transfer reaction rate constant can be calculated by using the detailed 
balance relation and can be expressed in the form kji = kij exp(— / ksT) . 

The quantum yields $m of electronic escape via branch L, M and the quantum yields 
<3>g of direct ground state recombination can be characterized for 5-sites sequential kinetic 
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model by the expressions 



2T 

Pg = -^Pi(s^0 + ), (17a) 
2F 

0l = -^A(s-O + ), (17b) 

or 



where the expression + $jv/ + $g = 1 have to by fulfilled. The analytical expressions for 
the ratio of the quantum yields have the forms 



<Pl fc 13 fc 35 (^f fc 24 + k 21 (k 42 + ))gfr 
0l _ k n k 35 ^ 



2Tl 

(18b) 



(^35 + M^3 + ^))^' 

The results of numerical calculations of QY's rate constants for the sequential model in 
both branches of RC for different samples of RC are collected in Table I. 

The expressions for the electron transfer are given by the inverse Laplace transformation. 
Therefore firstly we apply the Laplace transformation to P(t) in system of Egs. 16. Where 
the Laplace transformation is defined as 

POO 

P(s) = / e- st P(t)dt. (19) 
Jo 

Next we apply the inverse Laplace transformation of P(s) where the inverse Laplace trans- 
formation is represented by a set of simple poles of P(s). Evaluating it we obtain 

5 

P(t) = J2^ t } (20) 

3=1 

where aj are amplitudes and kj are rate kinetic constants describing the electron transfer. 

With using the model described above we would like to find kinetic of the reaction centers 
of Chloroflexus aurantiacus [17| where on the M-branch the BChljv/ is replaced by BPIim in 
corresponding position. Thus C. aurantiacus RCs contain altogether three BPh molecules 
and only one BChl monomer. To characterized C. aurantiacus we start from the set of 
parameters that characterize the kinetics of wild-type (WT) RCs of Rb.sphaeroides. We use 
the following values of input parameters: reorganization energies Ay = 800 cm -1 , electronic 
couplings Sij = 0.5 cm" 1 , high frequency modes lu^ = 1500 cmT x where i,j = 1,3,5 



for L-side and i,j = 1,2,4 for M-side of RC. The values for electronic couplings V 2 4 = 
V35 = 32 cm" 1 , V12 = V13 = 20 cm' 1 were used. The sink parameters 2T M /h = 2T L /h = 
(200ps)~ 1 ,2Ti/h = (170ps) _1 were used in accordance with experimental observation which 
characterize the ET to quinone molecules and decay to the ground state. Because of BChl^ 
is replaced by BPh^ in corresponding position we decrease the free energy in site 2. The 
calculated rate constants and quantum yields for the concrete energy levels are collected in 
Table I. We get the following occupation probabilities for different sites: 

P^t) = 0.05e- - 6t + 0.15 e - a37t + 0.8e- a22t , 

p 2 (t) = -0.25e~ 0M + 0.05e-°- m + 0.2e~ - 22t , 

P 3 (t) = -0.03e- a6 ' - 0.96e-°- 37 ' + 0.99e- a22t , 

P 4 (t) = 0.2e-° M - 0.07e-°- 37 ' - 0.44 e -°- 22 ' + 0.3 e - acm + O.Ole" - 009 ', 

P 5 (t) = 0.02e- 0M + 0.8e- a37 ' - 1.52 e - a22t - 0.4 e -°- 01 ' + Lie"™. (21a) 

The exponential components with very small amplitudes were neglected in the expressions 
above. Time evolution of the occupation probabilities is shown in the Fig.2. 




FIG. 2: The occupation probabilities P(t) for the reaction center Chloroflexus aurantiacus in the 
case if Vu = V13 . 

We can see that in this case we get electron transfer through the M branch, which is not 
in accordance with experimental observation. To avoid this discrepancy we must assume 
the asymmetry in the electronic coupling. To describe experimental kinetic of Chloroflexus 
aurantiacus RC we used the following asymmetry in electronic couplings: V12 = 10 cmT x 
and V13 = 15 cmr x . We weakly decrease the coupling constants in comparison with previous 
case, because of the kinetic in this RC is slower then in the WT RC. The calculated rate 

10 



constants and quantum yields for the concrete energy levels are shown in Table I, second 
line. We found 



Ps(t) 

p 4 (t) 

Ps(t) 



0.002e" a53t + 0.018e" ad4t + 0.98e" 
-0.05e- a53 * + 0.005e-°- 34 ' + 0.04 e -°- 12 ' + 0.005 e - aolt , 

-0.12t 



-0.34* 



-0.12* 



-0.4e-°- 34 ' + OAe 



0.053e 



-0.53* 



0.003e~ a34 * - 0.22e'°- 12t + 0.13e 



-o.oit 



0.04e 



-0.009* 



0.4e" 



-0.34* 



1.3e 



-0.12* 



0.5e 



-o.oi* 



+ 1.4e 



-0.009* 



In Figure 3 the behavior of the occupation probabilities Pi(t) is shown. 




20 40 60 80 100 



t[ps] 




80 100 



t [ps] 



(22a) 



FIG. 3: The occupation probabilities P(t) for the reaction center Chloroflexus aurantiacus in the 
case if Vi 2 / Vi 3 . 



Now we want to elucidate the electron transfer in H(M182)L. In this mutant BChl^/ 
is replaced with BPhi\/. The new cof actor is referred to as 6m- It is reasonable that in 

n 

the H(M182)L mutant the state P + <Pm is lower in energy than P + BChl~ M in WT [22]. To 
explain the electron transfer in this mutant we started from incoherent model. In this model 
we assume that the energy of P + 4>m is lower than the free energy of the state P + BPh~ h/I 



231 ] . The value of the free energies used to calculate the rate constant are listed in Table I. 
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The occupation probabilities Pi(t) in the case of H(M182)L mutant are found in the form 



Pi(t) = 0.33e- a38 * + 0.67e-°- 25 ', 

P 2 (t) = -0.02e-°- 52t - 0.07e~ a38 * - 0.23e-°- 25 * + 0.32e 
P 4 (t) 



-0.0006* 



= -1.4e- a38 * + 1.4e 



-0.25* 



. 04e -0.52 t _ . 036e -0.38 t _ . 04e -0.25 t + .036 e - a ° 006 ' , 



Lie" - 38 * - 1.602e-°- 25 * + 0.5 e - a01t + 0.002 e - a0006t , 



(23a) 



and the behavior of the occupation probabilities is shown in Fig. 4. It was assume that the 
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FIG. 4: The occupation probabilities Pi(t) for the mutant H(M182)L in the case of incoherent 
model of electron transfer. 



free energy of P + 4>m * s significantly below P + BPh~ M because of the electron transfer stops 
at 4>m- We get a very small probability to find electron on BPh M but the quantum yields 
through the branch M is substantial(Table I). 

Now we intend to elucidate the observed ET kinetics with partially coherent models. It 
means, that we assume that the reorganization energy for ET from state P + 4>~ M to state 
P + BPh^j is practically zero. The electron kinetic have to be described by the following 
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system of equations 



2T 

d t Pi(t) = -(—± + k 12 + k Vi )P 1 (t) 

+k 21 P 2 (t) + k 31 P 3 {t), (24a) 

d t P 2 (t) = -k 21 P 2 (t) - [ W 24 (t - r)P 2 (r)dr + k X2 P x {t) + / W 42 (t - r)P 4 (r)rfT<24b) 

Jo Jo 

d t P 3 (t) = -(k 35 + k 31 )P 3 (t) + fc 13 Pi(^ + k 53 P 5 (t), (24c) 
^P 4 (t) = ^P 4 (t)-J W 42 (t-r)P 4 (r)dr + J W 24 (t - r)P 2 (r)dr, (24d) 



2r 

d t P 5 (t) = -(-^ + k 53 )P 5 (t)+k 35 P 3 (t). 



(24e) 



In this case the memory function W 24 = W 42 can be expressed in the form: W 24 (t) = 



27r%^Re| exp - exp H£ \ £4; t j, where T 2 = for our kinetic model. We now 

use this partially coherent model of RC to describe H(M182)L mutation of RC. The results 
of our numerical computations are collected in Table I. We found the following expressions 
for occupation probabilities Pi(t): 



r , M+r , 2 1 



Pl(t) 

Ps(t) 

p 4 (t) 

Ps(t) 



0.6e- a42t + 0.38e- - 3t + 0.02e- a002 ', 

_ .3e-°- 42 ' - 0.201e-°- 3t + 0.001 e -° 009t + 0.5 e - a002t + 0.0004 e - a005t sin(3H), 
_l.4e-o.42t + L4e -o.3 t> 

-0.03e- a42t - 0.02e- a3t + 0.03 e - a009t + 0.02 e - a002t - 0.0006 e - a005t sin(31f), 



= e 



"0.42* _ 1Ale ~0M + _ 3e -0.01 t + . 01e -0.009t + O . le -0.002*_ 



(25a) 



The behavior of the occupation probabilities Pi(t) is shown in Fig.5. The calculated values 
of parameters for coherent model which describe the mutant H(M182)L RC's are collected 
in Table I. 
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FIG. 5: The occupation probabilities P(t) for the mutant H(M182)L in the case of coherent model 
of electron transfer. 
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£2 


C3 


£4 


£5 


IA12 


1/foi 


1/% 


1/*31 


1/A&4 


1/&42 




Sample 


K 


cm~ l 


cm -1 


cm~ l 


cm -1 


ps 


ps 


ps 


ps 


ps 


ps 


$G 


C.aurant. 


295 


-50 


-450 


-1000 


-2000 


12 


15 


6 


52 


2 


188 


0.02 0.31 0.67 


V12 = Viz 


200 










15 


21 


5 


138 


2 


1559 


0.02 0.25 0.73 


C.aurant. 


295 


-50 


-450 


-1000 


-2000 


47 


60 


11 


93 


2 


188 


0.05 0.16 0.79 


V 12 + Viz 


200 










58 


84 


10 


246 


2 


1559 


0.05 0.12 0.83 


H(M182)L 295 


-1600 


-450 


-1000 


-2000 


8 


10005 


6 


52 


37 


2.1 


0.02 0.41 0.57 


Incoherent 200 










9 


443897 


5 


138 


132 


1.8 


0.02 0.37 0.61 


H(M182)L 295 


-850 


-450 


-1000 


-2000 


5 


292 


6 


52 


4394 


4394 


0.03 0.12 0.85 


Coherent 


200 










4 


1867 


5 


138 


4394 


4394 


0.02 0.38 0.6 



TABLE I: The computed rate constant 1/fejj and quantum yields dependent on temperature for reaction 
centers and some mutants of RC's. The rate constants l/A^s = 3(4)ps for T = 295(200)7^ and 
1/&53 = 5166(264544)ps for T = 295(200)-?T are the same for all RC's and mutations describing in 
the Table I. 



IV. CONCLUSION 



We have dealing with electron transfer in the reaction center of Chloroflexus aurantiacus 
and Rhodobacter sphaeroides H(M182)L mutated RC. In spite of their structural similarity, 
the functionality is very different. H(M182)L mutant reveal the M brunch active in electron 



transfer. In the previous papers 
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25 



26. 



27 



were discussions about what is dominant 
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factor which causes the asymmetry in the electron transfer. At the beginning it was assumed 
that the asymmetry in the coupling parameters is dominant. The later the experimental 
work brought a doubt about dominance of electron coupling as a mechanisms which cause 



the asymmetry in ET trough branches 



22 



23]. Now, it is assume that the asymmetry 



in energetics also contribute to asymmetry of ET through M and L branches. We have 
showed that in the Chloroflexus aurantiacus RC we must have minimally 2:3 ratio of the 
electron transfer integrals for P*<pM ^ P + 4>ai an d for P^BChl^ <-> P + BChlJ j to explain 
the observed ET kinetics. In the case of Rhodobacter sphaeroides H(M182)L mutated RC 
we used two models to describe the ET. In incoherent model we must use a very low free 
energy of the -P + M state in comparison to the free energy of P + BPhJ 4 state to get a small 
probability to find electron on the BPhu molecule. Despite this assumption we get relatively 
strong outlet through the M branch. In the partially coherent model we assume that the ET 
between molecules 0m and BPfiM have coherent character, this means that ET is so fast 
that the bath does not have sufficient time to relax to the new thermal equilibrium before 
the particle moves away. The result is outlet trough the M branch which is in accordance 
with experimental observations. The problem of both coherent and incoherent models is 
that predict not enough decay to the basic state. The similar model ought to be used to 
characterized electron transfer in PSII (PSI) system, where dimer molecules of BChl are 
not so close as in the bacterial reaction centers. The temperature dependence of the kinetics 
and QYs was also computed. We can see that there are differences in dependence of QY on 
the temperature in both models which were assumed. 

The work was supported in part by VEGA grant 2/7056/27. of the Slovak Academy of 
Sciences, by the Science and Technology Assistance Agency under contract No. APVT-51- 
027904. 



[1] J. Deisenhofer et ai, Nature 318 (1985) 618. 

[2] J.P. Allen and J.C. Willams, FEBS Lett. 438 (1998) 5. 

[3] J. Amesz, J. Photochem. Photobiol. 30 (1995) 89. 

[4] H. Sakurai, N. Kusumoto, and K. Inoue, Photochem. Photobiol. 64 (1996) 5. 

[5] J.H. Golbeck, Proc. Natl. Acad. Sci. U.S.A. 90 (1993) 1642. 



15 



[6] W.D. Schubert et al, Nat. Struct. Biol. 3 (1996) 965. 

[7] Van Brederode et al, Biochemistry 36 (1997) 6855. 

[8] B. A. Heller, D. Holten, and C. Kirmaier, Science 269 (1995) 940. 

[9] Ch. Kirmaier and D. Holten, Proc. Natl. Acad. Sci. U.S.A. 87 (1990) 3552. 

[10] E. Takahashi and C. A. Wraight, Biochemistry (1992) 855. 

[11] V. A. Shuvalov and L. N. M. Duysens, Proc. Natl. Acad. Sci. U.S.A. 83 (1986) 1690. 

[12] J. N. Gehlen, M. Marchi, and D. Chandler, Science 263 (1994) 499. 

[13] M.G. Midler, G. Drews and A.R. Holzwarth, Chemical Physics Letters, 258 (1996) 194. 

[14] M. Marchi et al, J. Am. Chem. Soc. 115 (1993) 4178. 

[15] S. Tanaka and R. A. Marcus, J. Phys. Chem. B 101 (1997) 5031. 

[16] M. Bixon, J. Jortner, and M. E. Michel-Beyerle, Chem.Phys. 197 (1995) 389. 

[17] M.G. Midler, K. Griebenow and A.R. Holzwarth, Biochimica et Biophysica Acta, 1098 (1991) 
1. 

[18] M. Sparpaglione and S. Mukamel, J. Chem. Phys. 88 (1988) 3263. 

[19] R. Zwanzig, Physica 30 (1964) 1109. 

[20] N. Hashitsume, P. Shibata, and M. Shingu, J. Stat. Phys. 17 (1972) 253. 

[21] J. Jortner, J. Chem. Phys. 64 (1976) 4860. 

[22] E.Katilius,T.Turanchik,S.LIn,A.K.W.Taguchi,N.W.Woodbury, J.Phys.Chem. B 103 (1999) 
7386. 

[23] Ch. Kirmaier, Ch. He, and D. Holten, Biochemistry 40 (2001) 12132. 

[24] M. Pudlak, R. Pincak, Physical Review E 68 (2003) 061901. 

[25] M.E.Michel-Beyerle,M.Plato,J.Deisenhhofer,H.Michel,M.Bixon and 

J.Jortner,Biochim.Biophys.Acta 932 (1988)52. 

[26] R. Pincak, M. Pudlak, Physical Review E 64 (2001) 031906. 

[27] M.Plato,K.M6bius,M.E.Michel-Beyerle,M.Bixon and J.Jortner, J.Am.Chem.Soc. 110 (1988) 
7279. 

[28] M. Pudlak, R. Pincak, Chemical Physics Letters 342 (2001) 587. 



16 



